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20.  Abstract 


are  presented.  Specifically,  it  is  illustrated  how  the  detection  performance 
of  the  technique  is  established  by  the  total  energy  in  the  spectral  set  rather 
than  by  the  levels  of  individual  spectral  elements  in  the  set.  Furthermore,  it 
is  shown  how  minimum  variance  spectrum  analysis  actually  can  suppress  components 
in  such  a set.  Finally,  simulated  examples  are  presented  which  illustrate  the 
performance  of  the  technique. 
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A TECHNIQUE  FOR  COHERENT  DETECTION  AND  RELATING  OF  SPECTRAL  SIGNALS 


INTRODUCTION 


In  many  frequency  and  wavenumber  analysis  problems  the 
objective  is  to  detect  multiple  spectral  terms  in  simultaneous  fashion. 
Typically,  this  detection  objective  is  implemented  by  power  spectral 
estimation  followed  by  a search  routine  designed  to  relate  detected 
components  according  to  some  a priori  determined  structure.  For 
example,  spectral  terms  could  be  related  according  to  harmonic  family 
in  frequency  analysis  or  a plane  wave  arrival  structure  with  possible 
multipath  in  spatial  wavenumber  analysis.  Until  recently,  the  detection 
aspect  of  this  procedure  has  been  performed  as  if  the  spectral  terms 
being  detected  were  characterized  by  random  phase.  As  such,  phase  in- 
formation is  ignored  and  detections  are  performed  sequentially  on  the 
basis  of  estimated  power  spectrum.  More  recently  1»2.3  sequential 
adaptive  schemes  have  been  considered  which  exploit  stationary  phase 
estimators  which  yield  coherent  detection  algorithms  that  exhibit 
improved  detection  performance.  However,  very  little  consideration  has 
been  given  to  exploiting  potential  coherence  between  the  complex 
envelopes  of  multiple  spectral  components  with  the  objective  of  improv- 
ing both  detection  and  relation.  Intuitively,  one  feels  that  the  best 
detection  scheme  should  have  implicit  an  estimation  algorithm  for 
simultaneous  coherent  detection  and  component  relating.  In  the  follow- 
ing sections,  such  an  approach  to  coherent  spectral  component  relation 
which  occurs  simultaneously  with  detection  is  presented.  The  proposed 
scheme  is  based  on  a principal  component  analysis  of  the  random 
Fourier  transformed  data  vector  4 and  more  specifically  an  orthogonal 
vector  analysis  of  the  cross-frequency  correlation  (CFC)  matrix.5  This 
analysis  decomposes  the  CFC  matrix  into  its  eigenvectors  and  indicates 
that  the  frequency  information  on  harmonic  sets  with  coherent  envelopes 
appears  in  a single  eigenvector. 

THE  CROSS- FREQUENCY  CORRELATION  (CFC)  MATRIX 


Consider  the  frequency-to-frequency  Hermitian  correla- 
tion matrix  R with  term  in  the  i-th  row  and  j-th  column  given  by 
s 
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where  ( ) is  the  statistical  expectation  operator,  ( )* 

indicates  the  matrix  complex  conjugate  transpose  and  X(i/,  p)  is  the  p-th 
time  sample  of  a complex  envelope  function  for  a narrowband  random 
process  centered  at  the  discrete  frequency  v . Two  important  examples 
of  this  type  of  complex  envelope  time  series  are 

I output  time  series  from  a phase  rate 
normalized  discrete  frequency  analyzer  at 
frequency  v 6 

or  (3) 

output  time  series  from  a beamformer 
with  steering  at  spatial  frequency  (wave 
number)  v for  a particular  temporal  frequency. 

The  CFC  matrix  R is  illustrated  in  figure  1.  The  diagonal  of  R is  the 
discrete  frequency  power  spectrum  of  the  random  process  X(i>).  In  most 
frequency  analysis  based  detection  systems,  estimates  of  this  power 
spectrum  obtained  with  appropriate  temporal  averaging  and  spectral 
smoothing  is  the  only  information  used  to  detect  the  presence  of 
possibly  multiple  spectral  components  in  the  interval  [i/0,  j] . 

An  extraction  scheme  which  uses  only  the  diagonal  of  R is  commonly 
referred  to  as  being  incoherent  since  it  ignores  both  the  possibility 
of  nonrandom  (or  at  least  short  term  stable)  phase  rate  in  the  term 
Xfi'i*  p)  and  possible  coherence  between  the  "complex  envelopes" 

X(wi,  p)  and  X(vj,  p)  at  frequencies  Vi  and  i/j  respectively. 

In  terms  of  improving  detection  performance,  a logical 
question  to  ask  is,  "What  off-diagonal  terms  in  R might  exist  which 
constitute  useful  information?"  In  answering  this,  consider  as  an 
example  the  spectrum  analysis  problem  where  the  signal  of  interest  con- 
sists of  a harmonic  set  with  fundamental  frequency  with  harmonics  at 
{mvi)m=j.  Clearly, if  the  signal  component  of  the  phase  rate  normalized 
envelopes  at  each  frequency  in  the  set  have  stable  amplitude  and  phase, 
the  complex  envelopes  would  exhibit  correlation.  In  fact,  the  amplitude 
of  the  off-diagonal  term  rj  mi  in  the  i-th  row  would  be  proportional  to 
the  product  of  the  amplitudes  of  the  fundamental  and  the  m_th  harmonic. 
If  the  fundamental  at  has  completely  random  amplitude  and  phase 
envelope,  then  it  might  still  be  characteristic  of  the  harmonic 
envelopes  to  exhibit  random  amplitude  and  phase  variations  that  are 
proportional,  that  is,  correlated  with  those  of  the  fundamental.  Thus, 
while  there  might  be  no  justification  for  phase  detecting  any  element 
in  the  harmonic  set  individually  because  of  randomness  in  the  complex 
envelope  phase,  there  might  certainly  be  justification  for  exploiting 
finite  crosscorrelation  between  elements  of  a harmonic  set. 
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A somewhat  similar  argument  to  the  above  could  be  applied 
to  detection  of  signals  in  a Doppler  compensated  wavenumber  spectrum. 
However,  rather  than  a mechanism  as  structured  as  a harmonic  set  leading 
to  envelope  correlation,  multipath,  for  example,  might  well  provide 
spectral  components  at  different  arrival  angles  (wavenumbers)  which  have 
not  been  completely  "decorrelated"  by  propagation  through  the  medium.  As 
such,  the  complex  envelopes  would  lead  to  nonzero  cross- (spatial)  fre- 
quency terms  off  the  diagonal  in  the  corresponding  R matrix. 

Accepting  the  premise  that  there  might  be  useful  off- 
diagonal  information  in  the  cross-frequency  correlation  matrix  R leads 
to  the  determination  of  a systematic  search  routine  for  examining  the 
R-matrix.  Furthermore,  given  that  an  efficient  search  routine  can  be 
found,  it  must  be  determined  how  best  to  interpret  the  findings  of  such 
a search  in  terms  of  a detection  and  component  relate  objective.  In  the 
next  section  an  approach  to  the  decomposition  of  the  R matrix  which 
leads  to  a strategy  for  both  search  and  interpretation  is  developed. 


DECOMPOSITION  OF  THE  CROSS -FREQUENCY  CORRELATION  MATRIX 

In  the  following,  upper  case  bold  face  letters  are  used 
to  indicate  complex  matrices  and  underlined  matrices  represent  column 
matrices  (vectors). 

Let  the  deterministic  complex  N-dimensional  column 
vector JE  ^ describe  the  position  (frequency)  of  a spectral  component  in  a 
complex  N-dimensional  sample  space.  Furthermore,  let  the  complex  enve- 
lope, Sfc,  account  for  possible  random  gain  and  additive  random  phase  on 
each  nonzero  element  of  E.^.  Finally,  consider  an  observed  N-dimensional 
complex  random  data  vector  X consisting  of  K < N spectral  components  and 
additive  noise 


K 

X = X)  sklk  + N,  (4) 

k=l 


where  N is  an  N-dimensional  complex  sample  vector  of  zero  mean, 
identically  distributed  and  statistically  independent  noise.  In  par- 
ticular, the  element  in  the  i-th  row  and  j-th  column  of  the  noise  CFC 
matrix NN*  is 

' (5> 

= j 

j • (6) 


0O’  1 


0,  i ^ 


4 
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It  is  advantageous  to  introduce  the  following  matrix 
quantities:  An  N (row)  by  K (column)  matrix,  E,  wherein  the  spectral 
location  vector  E]<  is  the  k-tb  column. 


[*>**  "•  «k]  • 


(7) 


The  columns  of  E will  be  assumed  to  be  linearly  independent  but  not 
necessarily  orthogonal. 

A K by  K correlation  matrix,  P,  for  the  complex 
envelope  factor,  wherein  the  element  in  the  i_tb  row  and  j_tb  column  is 
given  by 


[P]..  = s.s? 

L J ij  1 j 

= P, 


ij 


(8) 

(9) 


Finally,  an  N by  N identity  matrix  1^  is  defined. 

With  the  above  the  CFC  matrix  can  be  written  as 


R = XX* 

= ere*  + o§,N. 


(10) 

(ID 


It  is  now  desirable  to  orthonormalize  the  spectral  location  vectors  in 
terms  of  a set  of  new  orthonormal  vectors  j where  if  0 is  an  N 

by  K matrix 


0 = [?i$2  • ’ ’ #K  • 


(12) 


then 


* 

* 


0*0  = lK. 


A K by  K upper  triangular  matrix 


B 


bllb12 

b22 

0 


blk 


kk 


(13) 


(14) 
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can  be  found  using  a Gram-Schmidt  procedure  where  0 can  be  expressed 
as 


0 


EB. 


(15) 


That  is,  each  spectral  location  vector  can  be  expressed  as  a linear 
combination  of  the  orthonormal  vectors  4>£.  Of  course,  the  orthonormal- 
ization procedure  presented  above  could  have  been  used  had  the  vectors 
Ek  not  been  linearly  independent.  In  this  case,  an  orthogonal  set 
j ’t’Jl } pJ=  i would  have  been  established  where  the  integer  L < K is  called 
the  dimensionality  of  the  signal  (sub-)space  defined  by  the  set  ^ Ej^ ^ k=  1 - 
For  the  duration  of  this  discussion,  however,  the  spectral  location 
vectors  will  be  assumed  linearly  independent,  since  this  ensures  that 
the  Gram-Schmidt  coefficient  matrix,  B,  is  a square,  nonsingular  matrix 
for  which  an  inverse,  B_1,  exists.  Using  (15),  equation  (11)  becomes 


R = 0 B" 1 PB" 1 *0  * + a2lN.  (16) 

The  Hermitian  matrix  B_1PB_1*  is  now  written  in  diagonal  form  as 


B-lRB-1*  = yly* 


(17) 


where  y is  a K by  K matrix  consisting  of  the  orthonormal  eigenvectors 
of  BPB-1*  and  1 is  a K by  K diagonal  matrix  consisting  of  the  eigen- 
values a2  >_  a|  >_  . . . 1 of  BPB"  1 * . It  is  worth  noting  that  (16) 
represents  a coupling  of  the  information  contained  in  the  spectral 
position  vectors  with  the  complex  envelope  factor  correlation  charac- 
teristics. 

Using  (17)  in  (16)  and  post-multiplying  by  0y  yields 


R0y  = 0y  1 + OqIjJ. 


(18) 


Observing  (18),  it  is  evident  that  the  matrix 0y  contains  K of  the  N 
eigenvectors  of  the  CFC  matrix  R as  columns  with  corresponding 

..  , , 'll.  oo 

positive,  real  eigenvalues  Xj 


i2  + 

K 


°1  + o2  > A2  = o2  ♦ 


°0  1 


a2  given  as  the  corresponding  diagonal  elements  of  the  matrix 


x = x 


°o'k- 


(19) 


The  remaining  N-K  eigenvalues  of  R are  identically  o2  corresponding 
to  a root  of  the  CFC  matrix  of  multiplicity  N-K.  The  elements  of  X 
in  (19)  have  been  assumed  to  be  arranged  in  decreasing  order.  Finally, 
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it  is  seen  that  the  N by  K matrix  of  eigenvectors  corresponding  to  the 
K largest  eigenvalues  of  R is  given  by 

M = (20) 

= EBy.  (21) 

Thus,  the  eigenvectors  (columns  ofM)  are  seen  to  be  linear  combinations 
of  the  spectral  position  vectors  where  the  combination  coefficients  are 
functions  of  the  inner  products  of  the  position  vectors £*  E j through 
B and  the  complex  envelope  factor  correlation  characteristics  through  y. 

Some  examples  of  the  results  obtained  in  this  section 
are  given  in  the  following  section.  Specific  attention  is  given  to  the 
case  of  two  (K  = 2)  spectral  components  because  more  complex  spectra  can 
be  considered  as  extensions  of  this  basic  signal  spectrum. 


ORTHOGONAL  DECOMPOSITION  OF  THE  CROSS- FREQUENCY  CORRELATION 
MATRIX  FOR  AN  M COMPONENT  SPECTRAL  SET 

For  simplicity,  let  the  observed  data  vector  consist  of 
two  spectral  components  plus  additive  uncorrelated  noise.  If  the  nota- 
tion ajj  = k | E . is  introduced,  then  a two-step  Gram-Schmidt  procedure 
to  orthogonalize^  E = [E  i E 2]  as  in  (15)  yields 


1 
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where 


A = ( ana22  - lai2  I 2)  1/2 


Therefore,  the  matrix 


C = B"  1PB_  1 ’ 


has  eigenvalues  given  by  the  roots  of  the  characteristic  equation  of 
B-ips-l*  t which  is 


o4  - CanPii  + a22p22  + 2 Re(pi2a*2))o2 


+ (ana22  - | a 1 2 | 2 ) (pnp22  - Pi2!2)  = 


As  a simplification,  let  the  following  definitions  he  introduced 


for  i=j  (intracomponent  spectral  gain) 

aii  = { 1/2 

J (G.G.)  a for  i/j  (intercomponent  spectral  gain) 


for  i=j  (intracomponent  spectral  power) 


(S.S^)  8 for  i/j  (intercomponent  spectral  power) 
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where  a shall  be  called  the  spectral  overlap  factor  (|a|<l)  and  6 shall 
be  called  the  envelope  correlation  factor  (|b|<1).  The  parameter  G is 
a scalar  gain  factor.  Using  (27)  and  (28)  in  (26)  gives 


Notice  the  duality  between  the  envelope  correlation  factor  8 and  spectral 
overlap  factor  a.  Two  special  cases  with  respect  to  (29)  are  of  interest 
in  that  they  give  the  extreme  values  of  and  a\  where,  in  particular, 
either  the  intercomponent  spectral  gain  factor,  a,  is  zero  due  to  the 
noncoincidence  in  frequency  of  elemental  spectral  components,  or  the 
intercomponent  spectral  power  is  zero  due  to  the  lack  of  coherence  be- 
tween spectral  component  envelopes.  Accordingly,  there  results: 


Noncoincident  spectral  components  (a  = 0 for  El  and  §2  orthogonal): 


• 4 -* 
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(note  Ej  = E2  given  a = 1) 


Ai  = G j S i + G2S2  + ctq 


Thus,  at  one  extreme,  the  occurrence  of  two  spectrally  disjoint 
components  in  the  same  eigenvalue  (eigenvector)  stems  from  correlated 
complex  envelope  factors  and  at  the  other  extreme,  the  complete  separa- 
tion of  disjoint  components  occurs  due  to  the  lack  of  envelope 
coherence.  The  significant  result  from  the  preceding  development 
is  that  by  orthogonal  decomposition  of  the  cross-frequency  correlation 
matrix  R,  simultaneous  extraction  of  envelope  coherent  spectral 
components  can  be  achieved,  because  a linear  combination  of  the  com- 
ponent position  vectors  forms  a single  eigenvector  of  R which,  in  turn, 
can  be  examined  for  its  spectral  content. 

To  further  emphasize  the  significance  of  (19)  and  (21) f 
consider  the  case  of  K spectral  components  with  orthogonal  position 
vectors,  i.e.. 


E*E . = G.6..  (6..  = dirac  delta),  (41) 

-i-j  1 13 


and  complex  envelope  power 


S.  = 
1 


(42) 


t h 

Let  the  K spectral  components  consist  of  L sets  where  the  l~z  set 
{icfl. } contains  components  for  which  the  magnitude  squared  coherence 


that  is,  envelope  coherence  exists 


( 1 for  i,je% 

| 0 for  iefi^  and  jefi^  (43) 

only  for  components  within  a set. 
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It  can  be  shown  from  (19)  and  (21)  that 


°0,MC 


' 1/2 
i i 


/ 2 2 CGiS±D 1/2M±  for  l < 


= oQ  for  N - L < l <_  N. 


Thus,  it  is  suggested  that  detecting  the  £ spectral  set  be  accomplish- 
ed by  comparing  an  estimated  to  Oq  and  relating  components  within 
the  £_th  spectral  set  be  accomplished  by  interrogating Mj;,  to  determine 
which  components  of  the  set  (§]<}£  , are  present  in  the  linear 
combinations. 

The  above  case  pertains  to  the  situation  where  each  of 
the  spectral  sets  is  disjoint,  i.e.,  Qjj,  = 0 with  l t k.  For  the 
case  12^  fl  flu  / 0,  which  includes  partially  coherent  complex  envelopes, 
leakage  between  sets  occurs.  However,  inclusion  of  a spectral  component 
within  a given  set  will  tend  to  be  in  proportion  to  the  amount  of 
coherence  which  exists  between  that  component  and  the  components  within 
the  set.  Returning  to  the  case  of  2 spectral  sets,  it  is  useful  to 
consider  the  behavior  of  a^,  for  |a|  and  | 8 | in  the  interval  (0,  1). 
This  consideration  assumes  a and  B to  be  real  valued.  (See  Cox7  for  a 
discussion  of  beamformer  mismatch  which  can  be  related  to  the  spatial 
interpretation  of  the  variable  a.)  N.L.  Owsley®  addresses  the  problem  of 
signal  suppression  in  noise  canceling  array  processors  due  to  signal 
and  noise  correlation.  This  type  of  suppression  can  be  related  to  the 

variable  B.  Figure  2 presents  plots  of  ( ) and  o|  ( ) 

for  various  values  of  GjSi  and  G2S2  versus  envelope  coherence  6 for  a 
specified  intercomponent  spectral  gain,  a.  The  transfer  of  energy  from 
one  eigenvalue  to  another  as  the  two  spectral  components  become  either 
spectrally  disjoint  or  envelope  incoherent,  as  the  case  may  be,  is 
readily  apparent.  Figure  3 gives  the  ratio  of  eigenvalues  o\l (dB) 
versus  B for  fixed  a and  three  values  of  G1S1/G2S2  (dB) . 


EXAMPLES  OF  SPECTRAL  SET  SEPARATION 


This  section  presents  several  examples  of  spectral  set 
extraction  by  the  method  of  orthogonal  decomposition  of  the  CFC  matrix 
which  was  presented  in  the  previous  sections.  The  presentation  is  in 
terms  of  an  N by  N,  2-dimensional  square  array  0^  with  eigenvalue  number 
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Figure  2.  Distribution  of  Component  Energy  Between 
Eigenvalues  as  a Function  of  Spectral  Set 
Overlap  (a)  and  Envelope  Correlation  (B) 
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(i)  and  eigenvector  element  number  (j)  as  the  independent  variables, 
i le  spectral  level  plotted  in  the  3-dimensional  format  of  figure  4, 
for  example,  is  the  value 

°ij  = Ailmij l2»  (46) 

where  m^j  is  thej-t^  element  in  the  N by  1 eigenvector corresponding 
to  the  eigenvalue  of  the  CFC  matrix  R.  Thus,  since  the  i-th  eigen- 
vector is  a linear  combination  of  the  spectral  position  vectors  for 
all  components  with  spectral  overlap  and/or  envelope  correlation,  the 
term  |mjj|2  will  give  an  indication  of  the  presence  of  energy  in  the  i_th 
spectral  set  at  the  j-th  frequency  cell . The  eigenvalue  weights  the 
normalized  component  energy  by  the  total  energy  in  the  harmonic  set. 

The  examples  which  follow  are  described  as  labeled  on  each  figure.  A 
CFC  matrix  for  N=64  frequencies  is  treated  in  each  case.  Figures  4 to 
8 present  cases  where  the  CFC  matrix  for  spectral  components  with  speci- 
fied phase  randomness  is  stipulated  analytically  and  formulated  for 
simulation  on  a general  purpose  computer.  Phase  randomness  is  specified 
in  terms  of  a uniform  distribution  of  phase  over  one  of  three  intervals: 


(0°,  0°)  coherent  (C) 

(0°,  180°)  partially  coherent  (PC) 

(0°,  360°)  incoherent  (I). 

When  a component  is  indicated  as  being  a member  of  a specific  harmonic 
set,  this  implies  that  the  phases  are  coherent  with  respect  to  other 
members  of  the  set,  that  is,  the  phase  difference  between  the  component 
envelopes  is  a constant  even  though  the  absolute  phase  of  each  component 
may  be  randomly  distributed. 

Figures  4 and  5 give  an  example  of  a three  set  harmonic 
extraction  operation.  All  parameters  for  the  two  runs  are  identical, 
except  for  a relative  decrease  of  6 dB  for  the  results  shown  in  figure 
5.  Notice  that  the  results  are  substantially  the  same  except  for  some 
artifacts  in  figure  5 due  to  operation  at  low  SNR,  i.e.,  small  eigen- 
value magnitudes.  Figure  6 has  the  same  parameters  as  the  case  of 
figure  4 except  components  2 and  6 are  partially  coherent.  This 
coherence  of  phase  leads  to  the  appearance  of  considerable  spectral 
energy  from  these  components  in  the  first  eigenvalue  and  corresponding 
vector.  Figure  7 illustrates  an  example  of  two  sets  with  components 
which  are  spectrally  coincident  but  incoherent  at  frequency  cell  30. 

This  results  in  the  appearance  of  energy  in  eigenvectors  1 and  3.  Also 
for  this  example,  the  partial  envelope  coherence  of  the  component  at 
frequency  45  results  in  a spreading  of  energy  over  several  eigenvectors. 
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Figure  5.  Spectral  Set  Separation 
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Figure  6.  Spectral  Set  Separation 
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essentially  as  predicted  by  the  results  of  the  section  on  ORTHOGONAL 
DECOMPOSITION  OF  THE  CROSS-FREQUENCY  CORRELATION  MATRIX  FOR  AN  M COM- 
PONENT SPECTRAL  SET.  Figure  8 contains  4 sets  with  spectral  coincidence 
at  frequency  cell  30  between  sets  1 and  4.  This  coincidence  leads  to 
coupling  between  spectral  position  vectors  as  predicted.  The  results 
illustrated  in  figure  9 derive  from  an  actual  random  frequency  "jitter" 
at  cells  15  and  45.  The  random  frequencies  are  stable  with  respect  to 
each  other  and  therefore  are  contained  in  the  same  eigenvectors.  This 
would  be  an  example  wherein  coherent  processing  of  each  component  in- 
dividually would  be  an  exercise  in  futility,  whereas  joint  processing 
in  terms  of  eigenvectors  would  yield  nearly  a 3 dB  increase  in  the  cor- 
responding eigenvalue. 

RELATION  OF  ORTHOGONAL  CFC  MATRIX  DECOMPOSITION  TO 
MINIMUM  VARIANCE  SPECTRUM  ANALYSIS 


A minimum  variance  estimate  of  the  complex  envelope  for 
a possible  component  of  the  random  data  vector  X with  spectral  position 
vector  D requires  the  minimum  variance  filter  for  X given  by 

H = G , (47) 

D*R1D 


where  D*D  = G.  The  expected  output  power  for  the  minimum  variance 
filter  is  given  by 


|H*X|  2 
G? 

D*R_1D 


(48) 

(49) 


In  the  preceding  sections  the  cross- frequency  correlation  matrix  R has 
been  decomposed  according  to 

R =M£M*  + oglN,  (50) 

which  can  be  written  as 


where 


J54k  = EByk  • (52) 

(Note:  A,  implies  the  column  and  A,  the  k'th  row,  respectively, 

of  the  matrix  A. ) 

The  inverse  of  R can  be  expressed  as 


M£(  1 + I k )-JM* 


M'PM* 


£ vkBSkMk 


k=l 


(53) 

(54) 

(55) 


where  f is  a K by  K diagonal  matrix  with  the  k-th  diagonal  element  given 
by 


¥ 


k 


(56) 


Thus,  the  output  power  of  the  minimum  variance  filter  is 


XXl9*Mkl2 

k=l  k K 


(57) 


Consider  the  case  of  two  spectral  components  with 
orthogonal  position  vectors,  i,e.,  a = 0,  (57)  becomes 


C - 


GSi  + ai 


^ ^mv  I 


l9*li|2 


s2 


GS2  + On 


|d*e2| 


incoherent 

envelopes 

(6  = 0) 

(58) 

coherent 


GSj  + GS2  + Oj 


|d*(s1/2E!  + s|/2E2) 1 2 envelopes 


(6  = 1) 
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An  important  special  case  of  the  above  exists  for 


ETE.  = G 

-i-i 


= **  ie 
= G p . | e^  i 

i . 1 i 


where  0 <_  | p ^ | <^|  1 where  | p.  | is  termed  the  generalized  cosine  of  the 
angle  between  the  vectorsEf  andE^.  Equations  (58)  and  (59)  can  now  be 
written  as 


(Si|N) 

1 + (SjN) 


(s2|n)  io 

1 + (S2 1 N)  Pz  = 0) 


- 1 . (S;  |n)  . (S2|N)  [<silN)l».l2  * (si  In) |p2!2 

+ 2(S1|N)1/2(S2|N)1/2|P1|  |p2|  COS  ( 0 1 - 02)  ] (6  = 1) 

where  (S./N)  = GS^/a^  is  the  postfilter  signal-to-noise  ratio  for 
8 = 1 in^he  absence  of  any  other  component  and  where  (Oj  -02)  is  the 
phase  difference  between  the  filter  outputs  D*Ej  and  D*E2.  Figures  9 
and  10  present  plots  of  the  ML  beamformer  output  power  for  two  spectral 
components,  where  the  two  components  exhibit  both  incoherent  envelopes 

(MV I , ) and  coherent  envelopes  (MVC,  ).  In  addition,  the 

conventional  filter  output 


|ycl2  - lQ*x| 2 , 


Ga§  C(Sj  |N)  |P1|2  + (S2 1 N)  | p 2 1 2 + 1)]  (6  = 0) 

Gag  [(Sj  |N)  |Pl  |2  + (S2  | N)  | p2 1 2 

+ 2(s1|N)1/2(s2|N)1/2|p1| |p2|  COS(01  - e2)  ♦ 1]  (8  = l) 

(64) 
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is  given.  In  figures  10  and  11  the  filter  output  power  is  plotted  as  a 
function  of  the  generalized  cosine  |pj|  of  the  angle  between  D*  and  E i 
and  | p 2 1 is  assumed  to  be  zero,  that  is,  D and  E?  are  considered  to  be 
orthogonal.  The  mechanism  for  suppression  of  signal  with  position 
vector  E j due  to  coherence  as  a function  of  (S2JN)  for  the  signal  with 
positioiTvector  E2  even  though  e 1 and  E 2 are  orthogonal  can 
observed  in  (62). 


In  the  above,  it  has  been  shown  that  for  incoherent 
spectral  component  envelopes,  minimum  variance  filtering  resolution 
sensitivity  for  the  position  vector  e 1 with  respect  to  the  desired 
filter  vector  D is  superior  to  conventional  filtering.  This  can  be 
easily  be  seen  by  considering  the  sensitivity  expressions 


8IW2 


9 I P 1 


1 + (Si|N) 


[1  + (1  - |piI2)(Si|n)]2 


Go2(Si|N)  (B  = 0)  (65) 


P2=0 


and 


9ly„ 


3 |p  1 


= Go2(Sj |N) , (6  = 0 and  B = 1)  (66) 


P 2 


=0 


for  | p 2 I = 0.  In  the  case  of  coherent  spectral  term  envelopes,  however, 

1 + (SjN)  + (S2|N) 


3I>W|2 


9 1 p 1 


[1  + (1  - |pi12)(Si|n)  ♦ (S2|N)] 


Go2  (SjN) 


2 0 


p2=0 


(B  = 1).  (6?) 


In  which  case,  as  the  signal-to-noise  ratio  of  the  second  spectral 
component  (S2/N)  increases,  the  sensitivity  to  S!  decreases  according  to 


’l^l2 

Ga2  (S ! |N)  |pi|2 

QL 

3 (S2/N) 

1 

II 

0 

CM 

a 

’l  + (1  - |pi|2)(Si|N)  4 (S  2 |N)]  ^ _ 

(68) 
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rr.v 


which  is  always  negative.  Thus,  it  is  apparent  that  the  application  of 
maximum  likelihood  filtering  to  coherent  envelope  spectral  components  in 
noise  can  give  poor  results.  Furthermore,  if  a spectrum  rich  in  coherent 
(and  partially  coherent)  components  is  to  be  analyzed,  a direct  decom- 
position of  the  CFC  into  orthogonal  components  could  provide  an  effective 
spectrum  analysis  technique. 


CONCLUSION 


A spectrum  analysis  technique  based  on  an  orthogonal 
decomposition  of  the  cross-frequency  correlation  matrix  has  been  pro- 
posed and  compared  in  some  respects  to  both  conventional  linear  and 
minimum  variance  spectrum  analysis.  Orthogonal  component  spectrum 
analysis  exhibits  the  useful  property  of  simultaneous  extraction  of  a 
set  of  spectrally  disjoint  components  which  are  characterized  by  co- 
herence between  the  envelopes  of  the  components  within  the  set.  It  has 
been  shown  that  minimum  variance  spectrum  analysis  for  this  class  of 
signals  can  actually  lead  to  worse  performance  than  conventional  analysis 
due  to  suppression  of  coherent  components.  Orthogonal  decomposition,  on 
the  other  hand,  exhibits  detection  performance  governed  by  the  ratio  of 
the  total  signal  power  in  an  envelope  coherent  spectral  set  to  the  noise 
level  in  a single  spectral  resolution  cell.  As  such,  orthogonal  spec- 
trum analysis  would  have  application  to  the  extraction  of  a limited,  but 
perhaps  interesting,  class  of  signals.  It  is  noted  that  the  method  of 
estimating  the  frequency  parameter  from  the  estimated  CFC  eigenvectors 
has  not  been  addressed  herein.  It  is  suggested  that  each  eigenvector 
could  be  interrogated  sequentially  using  a variety  of  techniques  such  as 
either  maximum  entropy  or  Pisarenko  spectral  decomposition  techniques. 
Reference  (9)  is  suggested  as  a survey  of  such  techniques.  Finally, 
adaptive  algorithms  based  on  gradient  search  techniques  are  available 
for  realizing  the  orthogonal  component  spectrum  analysis  technique  de- 
scribed herein.10  These  algorithms  will  circumvent  the  necessity  of 
first  estimating  and  storing  the  CFC  matrix  with  subsequent  orthogonal 
decomposition. 
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